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ABSTRACT 

We applied the Maximum Likelihood method, as an image reconstruction al- 
gorithm, to the BAT X-ray Survey (BXS). This method was specifically designed 
to preserve the full statistical information in the data and to avoid mosaicking 
of many exposures with different pointing directions, thus reducing systematic 
errors when co-adding images. We reconstructed, in the 14-170 keV energy band, 
the image of a 90x90 deg^ sky region, centered on (RA,DEC)=105°,-25°, which 
BAT surveyed with an exposure time of ~ 1 Ms (in Nov. 2005). The best sensitiv- 
ity in our image is ~ 0.85 mCrab or 2.0 x 10~^^ erg cm~^. We detect 49 hard X-ray 
sources above the 4.5 a level; of these, only 12 were previously known as hard 
X-ray sources (>15keV). Swift/XRT observations allowed us to firmly identify 
the counterparts fo r 15 objects, while 2 objects have Einstein IPC counterparts 



( iHarris et al.lll990l ): in addition to those, we found a likely counterpart for 13 ob- 



jects by correlating our sa mple with the ROSAT All- Sky Survey Bright Source 



Catalog (jVoges et al.lll999l ). 7 objects remain unidentified. Analysis of the noise 
properties of our image shows that ~ 75% of the area is surveyed to a fiux limit 
of ~lmCrab. This study shows that the coupling of the Maximum Likelihood 
method to the most sensitive, all-sky surveying, hard X-ray instrument, BAT, is 
able to probe for the first time the hard X-ray sky to the mCrab fiux level. The 
successful apphcation of this method to BAT demonstrates that it could also be 
applied with advantage to similar instruments like INTEGRAL-IBIS. 



Subject headings: galaxies: active - surveys - X-rays: binaries - X-rays: galaxies 
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Introduction 



More than 40 years after its discovery, the nature of the Cosmic X-ray Background 
(CXB) is still debated. Population synthesis models, based on unified AGN s chemes, explain 
the CXB spectrum us ing a mixture of obscured and unobscured AGN (e.g. 
1995l : lGiUi eTaPboOlh . 



Comastri et al. 



According to these models, most AGN spectra are heavily absorbed, and about 85% of 
the radiation produced b y super massive black hole accretion is obscured by dust and gas 
( jFabian fc Iwasawalll999l ). 



Deep soft X-ray surveys (0.5-2. keV) were able to resolve the majority (~ 80%) of 
the CXB flux into discrete sources (IHasinger et al.l Il998f). However the resolved fraction 
decreases with energy, being ~50-60% in the 6-8 keV band ( Giacconi et al.|[2002 : Rosati et al. 



20021 ) and even less above >8keV; the missing CXB component has a spect ral shape that 



is con sistent with a population of yet undetected, highly obscured AGN (see IWorsley et al. 
2005h . 



It is important to realize that highly obscured objects are detectable in X- rays only above 



10 ke V. Moreover, most of the energy of the CXB is emitted around 30 keV (iMarshall et al. 



19801 ) and the exact nature of the source population responsible for the background at these 
energies is unknown primarily because of the low sensitivity of previous X-ray telescopes 
operating above 15keV. 

All these reasons together motivate more sensitive observations of the hard X-ray sky. 



The B urst Alert Telescope (BAT; iBarthelmy et al. I l2005h . on board the Swift mission (IGehrels et al 
20041 ). launched by NASA on 2004 November 20, represents a major improvement in sensi- 
tivity for imaging of the hard X-ray sky. BAT is a coded mask telescope with a wide field of 
view (FOV, 120° x 90° partially coded) aperture sensitive in the 15-200 keV domain. BAT's 
main purpose is to locate Gamma- Ray Bursts (GRBs). While chasing new GRBs, BAT sur- 
veys the hard X-ray sky with an unprecedented sensitivity. Thanks to its wide FOV and its 
pointing strategy, BA T monitors continuously up to 80% of the sky every day. Early results 
from the BAT survey (IMarkwardt et al.ll2005l ) show that BAT is already ten times more sen - 
sitive than the previous hard X-ray all-sky survey performed by HEAO-1 (jhevine et al. 1984). 



Coded mask telescopes are, until the advent of next generation hard X-ray focusing op- 
tics, among the most sensitive instruments able to image the sky in the hard X-ray domain. 
Objects in the FOV cast part of the mask pattern onto the detector plane. Since the sources' 
signal is coded by the mask onto the plane this phase is also referred to as coding phase. 
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Thus, a decoding procedure is required in order to reconstruct the original sky image. A 
variety of methods can be use d to reconstruct the sky image in the case of a coded mask aper- 



ture (see iSkinner et aLl 119871 . for a general discussion on reconstruction methods). Among 
them, standard cross correlation of the shadowgram with a deconvolution array, the mask 
pattern, via FFT transforms, is the most often used. Generally, sky images are obtained for 
each individual observation, where an observation is defined as a period during which the 
attitude is stable and constant. Subsequently, another procedure, such as resampling and 
reprojecting, is needed in order to assemble the final all-sky image. 



Most of the extragalactic sources are very faint in the hard X-ray band. Thus their de- 
tection is challenging and requires sensitive techniques. We here describe the application of 
an alternative method which was designed to improve the sensitivity avoiding some of the 
disadvantages of the standard mask unfolding technique. 



This study has been performed in the framework of a campaign for optical spectroscopy 
analysis of a sample of "hard X-ray selected" extragalactic sources aimed at identifying new 
Sy2 galaxies. This paper discusses the meth od used to reco nstruct the survey image and 
presents the source catalog. A second paper (IRau et al.l 120071 ) describes in details the opti- 
cal campaign and the source identification proc ess; the spec t ral an alysis and the statistical 
properties of the source sample are discussed in lAjello et al.l (120071 ). 



The structure of the paper is as follows. In Section [2] we present details of the Maxi- 
mum Likelihood method that was developed to analyze the BAT data. In Section [31 we 
describe the analysis steps performed and we present and discuss the results of our image 
reconstruction algorithm. The last section summarizes the results. 



2. Spatial Model Fitting 



We apply "spatial model fitting", as described in [Strong et al.l (|2005[ ). to directly recon- 
struct the survey image from the raw detector data. "Spatial model fitting" means that a 
number of sky distributions, whose linear combination constitutes the model, are forward- 
folded through the full instrumental response in order to generate a model shadowgram. The 
model shadowgram, which is a linear combination of all model components, is then fitted in 
the full data space in order to get the most probable sky distribution. The actual search for 
unknown sources is then realized by moving a source probe in a grid over the sky. It is worth 
noting that no other steps, as image mosaicking, are required at the end of this process. 
This method was successfully applied to different kinds of experiments (i.e. COMPTEL and 
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INTEGRAL-SPI; IPiehl et al.lllQQSl : IStrong et all 120051 ) . Its development was driven by the 



capabilities of reducing systematic errors and noise e.g. the noise related with individual 
short images and the systematics when co-adding noisy images in the mosaicking procedure. 
This leads to an improvement in sensitivity over other methods, in particular reducing sys- 
tematic errors from background variations and resampling. The full information in the data 
is preserved and correctly treated in a statistical sense. 

The likelihood is the probability of the observed BAT data given the model. For our case it 
is defined as the product of the probability for each detector of each observation: 

^ = n Pi^''- 

ijk 

where 

Pijk = — ^ — • (2) 

is the Poisson probability of observing Uijk counts in pixel ij, during the A;-th observation, 
when the number of counts predicted by the model is %fc. 

The model is a linear combination of components; in the simplest case of 1 non variable 
source and 1 background component for each observation, we get: 

Oijk = Co X (A ® + Cfc X Bijk (3) 

where {A ® 5'"°''^'') is the convolution of the detector response (A) and a source of unit flux 
(S) at the sky position ao,6o and thus {A (g) S°'"'^°)ijk yields the prediction of counts from 
a unit flux source at the sky position ao, So in detector ij, during observation k; Btjk is the 
background prediction for pixel ij in observation k and cq and Cfc's are the parameters we 
want to estimate. 



For the analysis described in this paper, the background model comprises, for each 
observation, an empirical model (i.e. a 2 dimensional quadratic function similar to the one 
used by the tool batclean as described in section 13. ip and the model shadowgrams for all 
bright sources (see also section [3rT]) . The actual fit to the background is performed only once; 
during the source search only the normalization of the background in each pointing is allowed 
to vary. Since sources detected at this stage are faint, the background normalizations are 
expected to vary by very small quantities. Indeed, we verified that such variations were less 
than 10"'^ with respect to the background parameters determined before the source search. 

In the future, our method will allow to test more complex and physical background 
models (e.g. diffuse emissions). 
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2.1. Parameter model estimation 



The parameter values are found by maximizing the likelihood function, or, which is the 
equivalent, maximizing its logarithm: 



dlnL 



d 



, ijk ijk 



(4) 



where A is the vector of the parameters. This translates into the following set of equa- 
tions: 



dlnL sr^ , . I TLijk 



dlnL 

dck 



ijk 



ijk 



- 1 1 , VA; 



ijk 



(5) 
(6) 



which allows to estimate all parameters simultaneously. 

This set of equations can be solved only numerically and we use a modified Newton algorithm 
in order to find the solution. 



2.2. Source significance 

In the case of a single source component, the source significance can be estimated using 
the likelihood-ratio test. For this application, the null hypothesis is that no point source 
exists at the position under consideration and the background model can explain all the 
data. The alternative hypothesis is the converse. Two maximizations have to be done in 
order to calculate the likelihood Lq of the background (null hypothesis) and the likelihood 
of both source and background for the alternative hypothesis Li. The test statistic: 

T, = -2(lnLo-lnLi) (7) 



is expected, from Wilks's theorem (IWilksi Il938l ). to be asymptotically distributed as 
Xn ill the null hypothesis, where n is the additional number of free parameters that are 
optimized for the alternative hypothesis. Since in our case the source intensity is the only 
additional free parameter, the test statistics is expected to follow the xi distribution. Thus, 
the significance of a detection can be addressed as: 

S = ^Xlix)dx, (8) 
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which, after changing variables, become: 

Equation [9] is exactly the integral of the standard normal distribution from Ts to oo and 
so the significance of the detection is: 

T^/^a = ^/-2{\nLo-\nLl)a = na. (10) 

Hence, by definition, the significance fluctuations must be distributed as a normal Gaussian 
if everything is done correctly. 



2.3. Method implementation 

In case of large detector counts, the likelihood maximization is equivalent to the 
minimization, with the problem having the advantage that it can be solved faster analyt- 
ically. We have verified that in the case of large detector counts (>20) and large numbers 
of observations (>100) the two solutions are very similar and from now on we use the 
solution. 



The algorithm used is a parallelized implementation of spidiffit ( IStrong et al.l l2005l ) used 
for INTEGRAL-SPI data analysis. Parallelization was needed because of the size of the 
problem we are dealing with. The typical execution time needed to compute the analytical 
solution scales with n^ where n is the number of data points to fit (i.e.: number of BAT 
detectors, 32768, multiplied by the number of observations). A single minimization with 
2600 observations takes nearly 90s; the total execution time to generate a map of 450x450 
pixels would be ~200days. This time has been reduced to <15 days using an average of 15 
CPUs. We remark also that it is the first time that such an approach is applied to a problem 
of this large size. 



As shown in equation [3], the model is a linear combination of different components which 
can be specified at the input of the program. Source and background components are in 
general treated in different ways. Sources are assigned a single free parameter (their average 
intensity) while, as already discussed, the background components are allowed to vary from 
pointing to pointing. However, in case of variable sources, the user can specify that the 
source intensity is left as a free parameter in all pointings (or in time-contiguous groups of 
them). 
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We remark that for the analysis presented in the next sections, the program has been used 
in its simplest configuration, with only one constant source and normalizations of the per- 
pointing backgrounds allowed to vary. However, after the source search had been performed 
and source candidates identified we have used the ability to fit simultaneously all sources 
(each of which was again assumed constant in time). In fact, the simultaneous fit of all 
sources yields the best parameters (significances and fiuxes) and allows us to discard spu- 
rious detections. When the analysis is based on a large number of observations correlation 
( "cross-talk" ) between sources is negligible. 



2.4. Instrumental response 

As shown in equation [3l the first part of the model represents the source component 
(or components if more than one) and this is given in the most simple form by a point like 
source at position ao, 60 in the sky, forward- folded with the instrumental response. We have 
used a large set of Crab observations (> 1000) to develope a parametrized diagonal full 
instrumental response which enables us to predict the expected counts (essentially the term 
A (8> S"°'^" of equation [3]) from a unit fiux source as a function of energy and position in the 
field of view. The parametrized instrumental response was obtained in the following way 
(standard BAT software is reported in brackets): 



1. for each Crab observation a model shadowgram for the source position is computed 
(tool batmaskwtimg); 

2. for each Crab observation the source model and the standard background components 
of batclean (see Section 13. ip are fit to the data; 

3. the normalizations of the source components (in the different observations) are parametrized 

function of off-axis angle. 



The parametrized instrumental response is thus, for a given source position, the multipli- 
cation of the model shadowgram described in (1) and of a coefficient computed from the 
parametrization derived in (3). In this way the instrumental response accounts for the off- 
axi^ variation of the detected source intensity which the batmaskwtimg model does not take 
into account. The response, derived in this way, agrees with measured values to within 1 
sigma anywhere in the FoV. 



^The reader c an find more detail s about tlie off-axis variation of tlie source signal and other effects in the 
Appendix A.l of lAiello etaD |2007l ) 
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To improve the speed of the code during source search the full instrumental response 
was pre-computed over a 6' pitch grid in the whole BAT FOV. 

To conclude, in Fig. [1] we show the imaging reconstruction capabilities of our approach; 
two closeby faint sources, LMC X-1 and PSR B0540-69.3, are clearly detected in the image 
obtained using ~2600 observations. The good angular resolution of BAT is also preserved 
by our imaging reconstruction algorithm, in fact the two sources are separated by just 25' 
(for comparison the BAT Point Spread Function is 22'). 

3. Analysis 

In this section we describe the application of the Maximum Likelihood method to re- 
construct the image of ~ 1/8 of the sky using 8 months of BAT data. 

3.1. Data selection and screening 

We used 8 months of data, from April 2005 (when BAT data became public) to Novem- 
ber 2005. In order to secure optical follow-up with a dedicated observing campaign at La 
Silla, Chile in January 2006, we selected only observations with angular separation less than 
45 degrees from the zenith (RA=105 degrees, DEC=— 25 degrees). The all-sky analysis, still 
within the capabilities of modern super-computers, will be left to a future study. 

Swift-BAT survey data are in the form of 80 channels detector plane histograms (DPH) 
with typical exposure time of 300 s. 

In order to have a suitable clean dataset as input of the imaging reconstruction algorithm 
described in section [2], preprocessing must be carried out on the raw survey data. This 
preprocessing phase accomplishes two different goals: 1) data quality is monitored along 
the processing and 2) the very bright sources detected during each single observations are 
localized and inserted in the background model of the imaging reconstruction algorithm. 
The latter procedure can be justified as follows. The brightest sources (except the Crab 
Nebula) are known to be highly variable. However, since they are detected in general at high 
significance in a single observation, their intensities can be determined with good accuracy. 
Thus, inserting bright sources in the background model, rather than treating them as several 
independent components in the source model, allows to handle source variability in a natural 
way without increasing the size of the problem. 
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All the pre-processing was carried out using the latest available version of the Swift 
software contained in the HEASOFT 6.0.3. Below we report in brackets the name of the 
standard BAT tools used during our pre-processing. 

For each DPH, our pre-processing pipeline, does the following operations: 

1. data are rebinned in energy channels according to the gain-offset map generated on 
board [baterebin); 

2. the DPH is integrated along the energy axis, between 14 and 170 keV, and a detector 
plane image (DPI) is generated (batbinevt); 

3. a detector quality mask is created, where hot and cold pixels are masked out {bathot- 
pix). These pixels are identified as the wings of the distribution of counts for a given 
observation; in general 2% (and so roughly 1% on each side) of the distribution is 
excised; 



4. an empirical background model is fitted to the DPI (batclea 



5. the DPI and the background model are input to a FFT deconvolution algorithm which 
generates the sky image [batfftimage)] 

6. source detection takes place on the sky image and a catalog of all sources detected 
above S/N> 6cr is created (batcelldeteci); 

7. a model for each detected source is created and it is added to the background model 
of step 4. The source model is created using the measured source coordinates. These 
coordinates were preferred to the catalog position because of non-trivial systematic 
effects which produced a shift in the measured source coordinates as a function of po- 



sition in the FOV (see http://swift.gsfc.nasa.gov/docs/swift/analysis/bat_digest.html 
for more details) ; 

8. steps 4 to 7 are repeated until no new sources are detected in a single 300 s observation. 

In order to have the cleanest dataset possible we have applied cuts on the quality of the 
data. During the steps above data are screened on the basis of the following conditions: 



^the empirical background model built-in in batclean fits (for a given energy range) a quadratic spatial 
function plus a series of models which take care of detector edge effects for a total of 14 parameters. The 
user is also free to include sources or different background models. The reader can find more details about 
the batclean background model in the documentation included in the HEASOFT package (a copy of it is 
also available online at [http:/ /heasarc.nasa.gov/lheasoft/ftools/headas/batclean.html ) 
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• lock of the star tracker and pointing stability 

• spacecraft being outside of the South Atlantic Anomaly (SAA). This information is 
reported in the housekeeping data and is referred to a fiducial point inside the SAA. 

• BAT array rate <18000 counts s~^ 

• exposure being larger than 200 s 

• reduced of the background fit <1.5 

• > 9(7 detected sources must be within a distance of 0.1 deg from a known source 
otherwise they are thought to be spurious or transient. The observation is flagged for 
a later analysis, but not inserted into the final dataset. 

In table [1], we have listed the fraction of exposure which is rejected if a single data 
quality cut is applied to the data used for this analysis. For the current dataset, ~ 34% of 
the overall exposure time was rejected because the data did not meet one or more of the 
above mentioned criteria. 

After processing and screening the data according to these criteria, the final data set 
includes 2671 observations. These observations are input to the imaging reconstruction 
algorithm described in section [21 Fig. [2] shows the total exposure map of all pointings. 

All sources detected during the preprocessing phase are listed in table [2] along with their 
identification, their maximum and total significance (computed as the sum of the squared of 
significances) from this per-pointing analysis, and the number of detections. The distribution 
of the offsets of sources in table [2] from their catalog counterpart is reported in Fig. [31 The 
same graph shows the extremely good location accuracy of BAT which locates 95% of all 
sources, detected in single pointings, within 2'.2 radius. 

In order to understand the dependence of location accuracy on the source significance, we 
have analyzed all per-pointing detections, see Fig. [H and determined that the offset varies 
with significance accordingly to 

OFFSET = 4.94(±0.68) x (^^j^y^-^^^^-^^) [arcmin] (11) 
This analysis is based only on sources detected during individual pointings. 

3.2. Imaging reconstruction 

The 2671 DPIs along with their background models (created at step 4 in section [3TT1) are 
input of the imaging reconstruction algorithm. For this analysis we have used 1 parameter 
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for the source component; moreover we have allowed the normalization of the background 
component to vary separately in each pointing, leading to a total of 2672 parameters. The 
map is built in small segments of 5x5 deg^. A pixel size of 12' was chosen as the best 
compromise between computational time and the resampling factor of the PSF (~ 2 in this 
case). The significance image is shown in Fig. [51 



3.3. Setting the significance threshold 

There are several approaches in order to derive the best significance threshold. The 
Maximum Likelihood method leads to perfectly symmetric Gaussian, normal, noise in the 
pixels of the reconstructed image. Thus, the most straightforward approach is setting the 
threshold as the absolute value of the lowest negative fluctuation. In this case, since the 
negative fluctuations are given by noise, one should expect no false detection above this 
threshold. 

As it can be seen in the significance distribution reported in Fig. [6l no negative fluc- 
tuations larger than —4.3 are found. If we take into consideration the number of trials and 
the normal Gaussian distribution we get that above the threshold S/N ratio of 4.5 we expect 
a number of false detections of 0.7. We also made a Monte-Carlo simulation generating a 
large number (>1000) of sky images with Gaussian noise. We then counted all the excesses 
above the 4.5 cr level and found out that the number of expected false detection is 1.01, 
in agreement with the previous finding. A contamination of our sample of sources by ~1 
spurious detection was judged to be a good compromise between detection sensitivity and 
sample corruption (see Section for the chance connected to have a higher contamination). 
Hence we decided to fix the threshold to 4.5 a. 



3.4. Noise properties and Sky coverage 



The sky coverage is, for a given survey, the distribution of the survey's area as a func- 
tion of limiting flux. The knowledge about the sky coverage is particularly important when 
computing the number- flux relation (also known as LogN-L ogS distribution). We leave the 
derivation of the number- flux relation to a separate paper (lAjello et al.l 120071 ). but we are 
interested in deriving the sky coverage here as it brings crucial information about the sensi- 
tivity and noise properties of the survey. 
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The sky coverage as a function of the minimum detectable flux Fmm is defined as the sum 
of the area covered to fluxes fi < Fmm- 

N 
i 

where N is the number of image pixels and Ai is the area associated with each of them. 
We have followed two procedures to compute the sky coverage of our survey area: 

• the ML method produces a flux map and an error map as output of the fitting pro- 
cedure. In order to get the sky coverage we multiplied the error map by the 4.5 cr 
threshold S/N ratio and then counted the area as in Equation [T2l 

• we computed the local (flux) image variance using a sliding annular region whose 
internal and external radii were 5 and 30 pixels respectively. The noise of a given pixel 
is thus computed as the variance of the pixels contained in the annulus centered on it. 
The central pixels are excised so that the background does not include contamination 
from the source region. This map is a true representation of the noise in our image. 
Again, we multiplied this noise map by our detection threshold of 4.5 a and then 
counted the area as in Equation [121 

The sky coverage computed in both ways does not present any significant differences 
testifying that the error computed by the ML method is very close to (if not the same as) 
the real noise term of the sky image. In the left panel of Fig. [7] we report the sky coverage 
of the entire area and for the extragalactic portion of the sky (selected imposing | 6 |> 15°). 
As it can be seen from the sky coverage, >75% of the surveyed area is sensitive to fluxes 
~lmCrab and all of it to fluxes >2.0mCrab. The limiting sensitivity in our image is a bit 
less than 0.9mCrab (or 2.05x lO^^^erg cm~^s^^). 



The analysis of the pixel noise as a function of exposure time (reported in the right panel 
of Fig. [7]) shows that the survey sensitivity scales oc T^°'^ denoting that systematic errors 
do not dominate over statistical ones. We then compared our sur vey sensitivity to recent 



results from the BA T and INTEGRAL-ISGRI hard X-ray surveys flMarkwardt et all 12005 



Bassani et al.ll2006l ). In order to perform the comparison, we transformed the sensitivities 



provided by the authors in different bands to sensitivities in a common band (20-100 keV); 
the comparison, which is shown in table [3l is done in two ways: once taking into account 
the threshold S/N used by the authors in their work, and then also based on a common 
5 (T-equivalent sensitivity. The main result is that for 1 Ms of exposure, our survey is one of 
the most sensitive. 
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3.5. Source detections and fluxes 



Source detection on ttie reconstructed image is a straightforward process since signifi- 
cance and flux maps are direct results of our reconstruction algorithm. All not-neighbouring 
pixels which meet the criterion S/N > S / Nthreshoid are identified. However at this stage we 
have lowered our detection threshold to (an optimally chosen) 3.5 cr. Indeed, our procedure 
of using a pre-computed response over a 6' pitch grid on a 12' pixel-size map might produce 
a small loss in the reconstructed sources' fluxes and significances. In order to overcome this 
problem, we have generated, for all the above candidates, a 5' pitch grid map using the 
correct instrumental response (i.e. not pre-computed on a 6' pitch grid). One such map has 
already been shown in Fig. [H Only those candidates whose significance, as derived from 
the oversampled small maps, exceeds the 4.5 a threshold are kept in the sample and fit with 
the instrumental point spread function (the batcelldetect is used here) in order to determine 
the most accurate source parameters. This procedure allows us to recover the correct source 
significance and flux at the cost of a slightly larger number of false detections. Indeed, due 
to the increased number of trials the expected number of false detection is now 1.5. We 
remark that our map is one realization over many; thus there exists a not-zero probability 
that the number of false detection exceeds the (averaged) value estimated here. Our Monte 
Carlo simulation shows that the probability of getting a number of false detections of 2, 3, 
and 4 is respectively 0.21, 0.09, and 0.02. 



3.6. Detected sources 



We have detected 49 hard X-ray sources in our survey. Four of these sources are residu- 
als caused by imperfect modeling (and inclusion in the background model) of bright sources 
which are detected in individual DPHs. These 4 sources (LMC X-4, EXO 0748-676, Vel X-1 
and V* V1055 Ori) are still detected in the reconstructed image with a S/N of 20-40. 

In table HI we report the coordinates and fluxes of all 45 serendipitous objects detected 
above the 4.5 a detection threshold. 



We have correlate d our sources with the ROSAT All-Sky Surve y Bright Source Catalogue 
( Voges et al.lll999l ) in the same way as in Stephen et al. (2006). In Fig. [HI we report the 



number of BAT sources which have at least one ROSAT source within a given radius. Also, 
to understand the contribution of chance coincidences to these associations, we performed 
a Monte Carlo simulation using 5x10^ positions randomly distributed in our fleld. Due to 
non- uniformity in the distribution of ROSAT sources, the probability of a chance association 
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increases slightly towards negative Galactic latitudes. Taking into account the highest den- 
sity of ROSAT sources (for —40° <b< —20°), we get from Fig. [8] that using a radius of 300" 
for the identification of our sources will yield a probability of chance coincidence of 0.015 (1 
wrong identification overall). 

The same figure yields also information about the BAT point spread function location 
accuracy (PSLA), as the BAT uncertainty in the position dominates the ROSAT error. Thus, 
assuming that the ROSAT position is the "true" source position and considering only the 
ROSAT associations, we fitted an inverted Gaussian to the curve of Fig. M (see Fig. [9]); we 
derived that 95% and 99% of all spatial coincidences are within 3'.3 and 5' respectively. Thus, 
using a 5' radius for source identification yields the best compromise between probability of 
finding the BAT counterpart and chance coincidence. 

It is not surprising that <70% of our sample is correlated with the ROSAT catalog 
since photoelectric absorption might play an important role. Using the ROSAT catalog 
we achieved to identify 30 of our sources. These sources are ge nerally the brightest of ou r 
sample and they were already detected by previous observatories ( iMacomb fc Gehrelslll999l ). 



Using the same 300" error radius, we searched for spatial coincidences betwe en our sources 
and both the HEAO-1 c atalog of high en ergy sources (ILevine et al.l Il984l ) and the 2nd 
INTEGRAL-IBIS catalog (IBird et al.ll2006l ). We found that 2 sources were already detected 
in hard X-rays by HEAO-1 and 7 objects, including also the previous 2, by INTEGRAL. 
All these 7 sources were already detected at low energy by ROSAT. Two addit ional sources, 
3C 227 and V* BG CMi have an Einstein IPC counterpart jHarris et al.lll99oh. 3C 2 27 was 
also detected during a long (11 ks) ROSAT-PSPC observations (ICrawford fc Fabianlll995l ). 



Some of the new sources can be identified using the narrow field X-ray telescope (XRT) 
on board Swift. With its 5" position accuracy XRT is able to pinpoint the source counter- 
part in less than 2 ks. We requested and obtained 3 foUowup observations of our targets 
(J0732.5-1330, J0823.3-0456 and J0918.6+1617 ) and this allowed us to firmly identify the 
counterpart of those sources (lAjello et al.ll2006l ). Other sources (e.g. J0916. 4-6221, J0519.5- 
3240, J0505.8-2351 and J0 920.8-0805) were ob se rved by XRT as par t of the ongoing effort 
of the BAT all-sky survey jTueller et al.ll2005al lbl: I Kennea et al.ll2005l ). We also searched the 
Swift archive for XRT observations covering the fields of our sources. 

A total of 17 sources can be firmly identified thanks to XRT. The results of all the identifi- 
cation efforts using X-ray catalogs and XRT, are reported in table HI Details of all sources 
identified using XRT are given case-by-case in the next section. 

Using the sources with a known X-ray counterpart, we report, in Fig. [TOl the sources' offsets 
from their catalog position as a function of significance. A fit to the data shows that the 
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offset varies with S/N as: 

OFFSET = 10.7(±1.9) x s/N-°-^^(^°-°^^ [arcmin] (13) 

Moreover, from the same plot we expect that for a 4.5 a detection the maximum offset be 
5'; this is in perfect agreement with what is shown in Fig. [91 

The offset derived for a 10 a source from the previous relation and from equation [11] (i.e. the 
same 10 a source is detected in the single 300 s sky image) is 1'.18 and 1'.26 respectively. The 
small difference between the per-pointing location accuracy and the accuracy in the summed 
image is due to the fact that Equations [11] and [13] are computed for different ranges of S/N. 
Indeed, sources detected in the survey image (sum of 2671 shorter observations) span the 
4.5-10 range of significance while most of the sources detected in single pointings have S/N 
greater than 9 a (see Fig. [1]). Thus, we can affirm that our survey analysis preserves the 
good location accuracy of BAT. 



3.7. XRT observations 

SWIFT J0407.6+0336. XRT observed this source field for 7ks on Jul 11, 2006. The 
only detected object, RA(2000)= 04 07 16.2 Dec(2000)=+03 42 24.3, is coincident with the 
Sy2 galaxy 3C 105.0 and distant 4'.7 from the BAT position. It is the first time that 3C 105.0 
is detected in X-rays. 

SWIFT J0505.8-2351. XRT observed this source field for 3.2 ks on Aug. 20, 2005. Only 
one source is detected within the BAT error box at RA(2000)=05 05 45.4 Dec(2000)=-23 51 
16.8 coincident with the Sy2 galaxy 2MASX J05054575-2 351139. This object was already 



identified as the BAT counterpart in iTueller et al.l (l2005bl ) . 



SWIFT J0519.5-3240. A 7ks XRT observation was performed on Nov. 26, 2005. In 
the XRT field only two objects are detected. The brighter one at RA(2000)=05 19 35.5 
Dec(2000)=-32 39 22.4 is only I'.l from the BAT position. The fainter one is detected at 
RA(2000)=05 19 25.8 Dec(2000)=-32 42 32.3 and it is only a 2a detection. The bright source 
is associated with the nearby Syl.5 galaxy ESO 362- G 018 (also detected by ROSAT as 
IR XS J051936. 1 - 323910 ). ESO 362- G 018 was already identified as the BAT counterpart 



bv lTueller et al.l fl2005bf ). 



SWIFT J0522. 6-3625. XRT observed this field for 899 s in May 26 2005. Only one 
source is detected at RA(2000)=05 22 57.8 Dec(2000)=-36 27 29.7 at 4'.6 from the BAT po- 
sition. The XRT source is coincident with ESO 362-G021 a BL Lac object already detected 
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by R OSAT and XMM at lower energies and by BeppoSAX in hard X-rays ( iDonato et al. 
2005h . 



SWIFT J0539.9-2842. An XRT observation of 14 ks took place on Dec. 8, 2005. A 
faint source is detected at RA(2000)=05 39 09.3 Dec(2000)=-28 41 01.5 coincident with the 
z=3.1 QSO PKS 0537-2843 and distant 2^.5 from the BAT pos ition. The QSO was discovered 
in X-rays by the Einstein observatory (jZamorani et al.lll98ll ) and then studied by ROSAT, 
ASCA and lately by XMM. The BAT detection in hard X-rays is the first to date. 



SWIFT J0550.8-3215. A 9ks XRT observation took place on May 21, 2005. A very 
bright source is detected at RA(2000)=05 50 40.4 Dec(2000)=-32 16 15.5 distant 2'4 from 
the BAT position. The XRT sourc e is associated with a well known blazar PKS 0548-322 
already detected in hard X-ray (see lDonato et al.ll2005l ). The blazar is then the BAT coun- 
terpart. 



SWIFT J0552. 1-0727. During 9ks of exposure on Apr. 8, 2006, XRT detects a bright 
source located at RA(2000)=05 52 11.5 Dec(2000)=-07 27 24.2. The object is coincident 
with the well known Sy2 galaxy NGC 2110 and its position is only 0'.4 aw ay from the BAT 



detec tion. The detection of NGC 2110 in the 3-20 keV band by RXTE flRevnivtsev et al. 



20041 ) and the presence of no other source in the XRT field secure the identification of NGC 
2110 as the BAT counterpart. 



SWIFT J0640.0-2553. During 2.8 ks of observation on Mar. 23, 2006, XRT detects only 
one bright source at RA(2000)=06 40 11.8 Dec(2000)=-25 53 41.5 coincident with the Syl 
galaxy ESQ 490- G26 (already detected in soft X-rays by ROSAT as RX J0640 11-25536). 
The source position is 2'.5 distant from the BAT detection. 



SWIFT J0732. 5-1330, aka SWIFT J0732. 5-1331. XRT observed this field for 4ks 
on Apr. 28, 2006. Only one source is detected in the XRT field at RA(2000)=07 32 
37.7 Dec(2000)=-13 31 08.6. This source is coincident with an USNO Bl star J073237.64 - 
133109.0. The source was already identified as the BAT counterpart by lAjello et al.l (120061 ). 
Follow -up measurements in the optical determined that this source is a new intermediate 
polar (jWheatley et al.ll2006l . and references therein). 



SWIFT J0759. 9-3844. XRT observed the field of this source for 7ks. Three sources are 
clearly detected. The brightest of them is located at RA(2000)=07 59 41.2 Dec(2000)=-38 
43 57.9 being only 0'.5 away from the BAT position while the remaining two are distant more 
than 10'. The brightest object is coincident with the INTEGRAL source ICR J07597-3842 
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and with the ROSAT source RX J075942.0 -384359. The fact that the o nly source within 



4' fro m the BAT position is also detected by INTEGRAL in hard X-rays (Iden Hartog et al. 



2004J ) makes this source the BAT counterpart. The probable association with an IR point 
source IRAS 07579-3835, and a 1.4-GIIz radio counterpart (NVSS archive) makes the case 
for a AGN nature of the object. This source was identified as being a Syl.2 galaxy during 
a recent optical spectroscopy foUowup (IMasetti et al.ll2006bl ). 



SWIFT J0823.3-0456. Only a single faint source is detected by XRT during 1.2 ks of 
exposure on Jan. 6, 2006. The source is located at RA(2000)=08 23 01 Dec(2000)=-04 56 
02.5 and 2' away from the BAT detection. The object is coinc ident with the gala xy FAIRALL 
0272 and was already identified as the BAT count erpart by lAiello et al.l (120061 ). An optical 
follow-up showed that the source is a Sy2 galaxy (IMasetti et al.ll2006al ). 



SWIFT J0918.6+1617. SWIFT J 0918.5+1617, aka SWIFT J0918.5+1618, is another 
source found thanks to our algorithm (jAjello et al.ll2006l ). During an XRT followup of 0.6 ks, 
the only detected source is located at RA(2000)=09 18 25.8 Dec(2000)=+16 18 20.8 (2'.5 
away from the BAT position) and coinci dent with the galaxy Mrk 704. Mrk 704 was previ- 



ously detected in soft X-rays by R OSAT (ISchwope et a. 



the galaxy was found to be a Syl (IMasetti et al.ll2006al ) 



2000). In a recent optical followup. 



SWIFT J0920.8-0805. An XRT observation of 8.5 ks took place on Dec. 10, 2005. Only 
one source is detected in the entire field. Its position, RA(2000)=09 20 46.0 Dec(2000)=-08 
03 21.8, is coincident with the Sy2 galaxy MCG-01-24-012 and dista nt 2^.2 from the BAT 
position. This object was already identified as the BAT counterpart by lKennea et al.l (120051 ). 



SWIFT J0945. 9-1421. An XRT observation of 11 ks took place on Jul. 8, 2006. The only 
source detected inside the BAT error box is located at RA(2000)=09 45 42.0 Dec(2000)=-14 
19 33.7. The source is coincident with the Syl. 9 galaxy NGC 2992 and is distant 1'.5 from 
the BAT detection. The source was already detected in soft X-rays by ROSAT as IRXS 
J094541.9-141927. 



SWIFT J0947.6-3056. XRT observed the source field for 10 ks on Dec. 9, 2005. Only one 
bright source is detected at RA(2000)=09 47 39.8 Dec(2000)=-30 56 55.4 coincident with 
the Sy2 galaxy ESO 434- G 040 and distant onl y 0^.4 from the B AT position. The galaxy 
was also detected in hard X-ray by INTEGRAL (IBird et al.l 120061 ). 
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4. Conclusions 



We have presented an application of the Maximum Likelihood method as a deconvolu- 
tion technique used to reconstruct the sky image when dealing with a coded-mask instrument 
like BAT. The main difference with other image reconstruction algorithms, such as the stan- 
dard cross- correlation technique, is that a sky distribution model is forward-folded through 
the full instrumental response and fit to the detector plane counts in order to derive the most 
probable sky image. This is realized in a single step including data from many pointings and 
thus no image mosaicking is required. This study was motivated principally by the capabil- 
ities of ML to: 1) preserve the full statistical information in the data and 2) to reduce the 
systematic errors connected to mosaicking techniques which other methods cannot avoid. 
This leads to an improvement in sensitivity over other methods. 



Moreover, this study is motivated by the need to use sensitive imaging techniques for the 
study of the hard X-ray sky. Although deep soft X-ray surve ys (0.5-2.0 keV) were able to 
resolve the majority of the CXB emission into discre te sources faasinger erZlllQQsl ). only a 
minor fraction of the CXB above 8keV is resolved (IWorslev et al.l 120051). F urthermore, the 
bulk of the CXB radiation is emitted around 30keV (IMarshall et al.l Il980l ) and the exact 
nature of the source population responsible for the background at these energies is unknown 
because of the low sensitivity of previous hard X-ray telescopes. The BAT coded mask de- 
tector, on board the Swift mission, represents a major improvement in sensitivity for imaging 
of the hard X-ray sky; thus, we tested our ML imaging algorithm on BAT survey data. This 
study was also complemented by an optical spe ctroscopy campa ign aimed at identifying 
BAT-discovered extragalactic hard X-ray objects (IRau et al.l 120071 ). 



The results presented in the previous sections can be summarized as follows: after screening 
our dataset for bad data as discussed in Section [3?T| the final survey image obtained using the 
ML method presents a perfect Gaussian normal noise. We detected 49 hard X-ray sources 
above the 4.5 a detection threshold. Only 12 were previously known as hard X-ray emitters 
(previously detected by INTEGRAL or HEAO-1). 37 are new sources detected by BAT due 
to our image reconstruction method. 

The correlation of BAT sources with the ROSAT catalog shows the extremely good loca- 
tion accuracy of the BAT instrument which is also preserved by our algorithm. Also it is 
worth noticing that ~30% of our sources are not correlated with the ROSAT objects; this 
is most probably due to the presence of photoelectric absorption in some of the new BAT 
sources. The analysis of the limiting flux as a function of pixel exposure (see Fig. [7]) for the 
reconstructed image sum of all observations, shows that systematic errors do not dominate 
over statistical ones and that BAT should be able to achieve, in the future, a sensitivity 
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of O.SmCrab with 3 Ms of exposure (if systematics remain at this level). The sky coverage 
shows that 75% of the survey is covered to flux ~1 mCrab and all of it to fluxes > 2.0 mCrab. 
All of this makes this analysis one of the most sensitive surveys ever performed in the hard 
X-ray domain. 



The optical spectroscopy identification o f the new sources and a discussion about the optical 
properties are left to a separate pap er jRau et al. 2007 ) while the statistical and spectral 
X-ray properties will be discussed in lAjello et al.l (120071 ). 
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Table 1. Fraction of rejected data due to any single criterion 



TYPE 


FRACTION 


Star tracker lock 


9.7% 


Pointing stability 


17.4% 


Outside SAA 


10% 


Exposure > 200 s 


7% 


Others 


<1 % 


All conditions 


34 % 
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Table 2. Sources detected during data screening 
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16.1 


40.7 


19 


LXB 


H 1608-522 


244.9860 


-15.6526 


25.5 


76.0 


24 


LXB 


Sco X-1 


247.9922 


-48.8173 


8.8 


11.5 


2 


HXB 


IGR J16318-4848 


250.2720 


-53.7497 


8.9 


19.6 


6 


LXB 


H 1636-536 


251.4481 


-45.6088 


13.1 


21.8 


6 


LXB 


GX 340+0 


255.1860 


-41.6435 


27.0 


38.9 


8 


HXB 


EXO 1657-419 


255.9813 


-37.8334 


50.2 


98.3 


22 


HXB 


V* V884 Sco 


256.4394 


-36.4345 


12.1 


13.2 


2 


LXB 


Sco X-2 



'^Note that the significance of the detection depends on the source intensity, exposure and 
on the position in the FOV. The "TOTAL S/N" was computed as the sum of the squared 
significances in each observation. 
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Table 3. Sensitivities of different hard X-ray surveys. 



Instrument 


Ref. 


Energy 




Sensitivity^ 


Equivalent 5cr Sens'' 






[keV] 




(@lMs [mCrab]) 


(@lMs [mCrab]) 


INTEGRAL 


Bassani et al., 2006 


20-100 


5.0 


0.8 


0.8 


BAT 


Markwardt et al, 2005 


15-195 


5.5 


1.3 


1.2^ 


BAT 


this work 


14-170 


4.5 


0.86 


0.95 



^The sensitivity is computed by considering the noise-exposure relation provided by the authors 
[e.g. right panel of Fig. [7]) evaluated at the threshold S/N they used to detect sources. 

'^Sensitivities are computed assuming a threshold S/N of 5 a for all instruments/surveys. 



'^The sensitivity reported by iMarkwardt et al.l (120051 ) is referred to the all-sky analysis. Their 
best 5cr sensitivity in high-latitude fields is ~1.0mCrab for IMs of exposure. 



Table 4. Detected hard X-ray sources 



SWIFT NAME R.A. DEC Flux S/N Exposure ID Oflfsct Type Catalogs* XRT position (ref.) 





(J2000) 


(J2000) 


(10-1-'- cest 




(100 ks) 




(arcinin) 








J0407.6+0336 


61.9178 


3.6517 


4.21 


5.2 


330 


3C 105.0 


4.7 


Sy2 


X 


04 07 16.2, -1-03 42 23 


J0426.4-5710 


66.6021 


-57.1775 


2.02 


4.5 


803 


lAXG J042556-571lt' 


3.5 


Syl 


R 




J0433.1+0520 


68.2982 


5.3374 


7.0 


10.4 


316 


3C 120'= 


1 


Syl 


R 




J0451.6-0349 


72.9205 


-3.8240 


3.73 


4.9 


416 


MCG -01-13-025'' 


0.9 


Syl.2 


R 




J0505.8-2351 


76.4674 


-23.8666 


2.62 


7.2 


816 


SWIFT J0505.7-2348 


1.7 


Sy2 


X 


05 05 45.4 -23 51 17.0 (1) 


J0510.8+1631 


77.7224 


16.5265 


7.37 


4.7 


231 


CSV 6150^ 


2.5 


Syl.5 


R 




J0514.0-4003 


78.5146 


-40.0558 


5.02 


10.4 


934 


4U 0513-40 


0.8 


LXB 


R 




J0516.0-0007 


79.0096 


-0.1332 


7.6 


6.1 


373 


QSO B0513-002'' 


2.5 


QSO-Syl 


R 




J0517.1+1633 


79.2839 


16.5605 


3.80 


4.5 


233 












J0519.5-3240 


79.8844 


-32.6720 


4.78 


8.2 


909 


ESQ 362- G 018'' 


1.1 


Syl.5 


R, X 


05 19 35.5, -32 39 22 (1) 


J0519.7-4545 


79.9460 


-45.7557 


4.48 


5.3 


904 


Pictor A*" 


1.5 


Syl 


R 




J0522.6-3625 


80.6581 


-36.4233 


3.65 


4.7 


927 


ESQ 362-G021'' 


4.6 


BLLAC 


R, X 


05 22 57.9, -36 27 29 


J0529.4-3247 


82.3541 


-32.7965 


5.15 


10.9 


904 


V* TV Colb 


1.3 


CV-DQ* 


R 




J0534.5-5801 


83.6470 


-58.0200 


1.8 


5.6 


800 


V* TW Pic'' 


2.1 


CV 


R, I 




J0539.0-6406 


84.7717 


-64.1148 


2.10 


6.6 


726 


LMC X-3'' 


2.1 


HXB 


R 




J0539.5-6943 


84.8917 


-69.7210 


4.48 


8.4 


681 


LMC 


1.4 


HXB 


R, I 




J0539.9-6919 


84.9878 


-69.3230 


6.30 


8.1 


681 


PSR B0540-69.3'' 


1.4 


Pulsar 


R, I 




J0539.9-2842 


84.9953 


-28.7029 


5.68 


4.6 


858 


PKS 0537-286'' 


2.5 


Blazar 


R, X 


05 39 54.1, -28 39 54 


J0550.8-3215 


87.7165 


-32.2610 


3.40 


6.1 


883 


PKS 0548-322'' 


2.4 


BL Lac 


R, X 


05 50 40.4, -32 16 15 


J0552. 1-0727 


88.0411 


-7.4554 


32.04 


33.7 


477 


NGC 2110 


0.4 


Sy2 


X 


05 52 11.5 -07 27 24 


J0558.0-3822 


89.5237 


-38.3799 


3.62 


6.7 


913 


LEDA 75476'' 


2.8 


Syl 


R 




J0640.0-2553 


100.0031 


-25.8931 


2.54 


5.3 


793 


ESQ 490- G 26'' 


2.5 


Syl.2 


R, X 


06 40 11.8, 25 53 41 


J0727. 5-2406 


111.8951 


-24.1039 


1.9 


5.3 


850 








R« 




J0728.6-2604 


112.1626 


-26.0696 


1.87 


4.7 


851 


V* V441 Pup'' 


3.5 


HXB 


R 




J0731. 5-1-0955 


112.8752 


9.9214 


3.30 


5.9 


578 


V* BG CMi 


1.3 


CV 


E 




J0732.5-1330 


113.1328 


-13.5037 


1.85 


5.9 


820 


SWIFT J0732.5-1331 


1.1 


CV 


X 


07 32 37.8, -13 31 07 (3) 


J0739.6-3144 


114.9127 


-31.7496 


2.77 


6.2 


846 












J0743.0-2543 


115.7501 


-25.7314 


1.86 


4.95 


869 








R° 




J0759.9-3844 


119.9822 


-38.7422 


4.90 


6.0 


804 


IGR J07597-3842 


0.5 


Syl.2 


R, I, X 


07 59 41.2, -38 43 57 


J0804.2-F0507 


121.0552 


5.1203 


4.1 


8.4 


736 


UGC 4203'' 


1.9 


Sy2 


R, X 


08 04 05.4 +05 06 49 


J0811. 5-1-0937 


122.8750 


9.6214 


2.26 


4.6 


700 








Rd 




J0823.3-0456 


125.8271 


-4.9401 


3.23 


8.1 


846 


SWIFT J0823.4-0457 


2.0 


Sy2 


X 


08 23 01.0, -04 56 02 (3) 


J0835.3-4510 


128.8308 


-45.1771 


17.1 


22.7 


705 


Vela PSR'' 


0.4 


Pulsar 


R, I, H 





C5 



Table 4 — Continued 



OVVi-T 1 iN/ViVlllj 


R A 




T71l1-ir 

r iUX 


C; /"NT 


Exposure 




wnsct 


lypG 


a 

Oatalogs 


jVixi positioii yxGx.j 






l^JZUUU ) 


(iO cgsj 








(arciiiin) 








J0838. 4-3559 


129.6151 


-35.9976 


4.01 


5.5 


791 


FRL 1146 


0.7 


Syl 


R, I 




J0839. 8-1214 


129.9556 


-12.2467 


2.43 


5.8 


866 


3C 206*^ 


0.4 


QSO 


R 




J0844.9-3531 


131.2411 


-35.5313 


2.97 


5.0 


780 












J0854.7-I-1502 


133.6828 


15.0371 


5.58 


5.1 


646 












J0917.2-6221 


139.112 


-62.359 


2.07 


4.5 


830 


LEDA 90443 


2.9 


Syl 


X 


09 16 08.9, -62 19 29.6 


J0918. 6-1-1617 


139.6505 


16.2987 


3.3 


6.7 


610 


MRK 0704 


2.5 


Syl.5 


R, X 


09 18 25.9, +16 18 20(3) 


J0920.3-5512 


140.0753 


-55.2135 


9.57 


15.6 


612 


4U 0919-54'' 


1.2 


LXB 


R, I, H 




J0920.8-0805 


140.2134 


-8.0872 


3.03 


6.4 


742 


MCG -01-24-012 


2.2 


Sy2 


X 


09 20 46.0, -08 03 21 (2) 


J0945. 9-1421 


146.4060 


-14.3007 


5.09 


5.5 


621 


NGC 2992*= 


1.5 


Syl.9 


R, X 


09 45 42.1, -14 19 34 


J0947.6-3056 


146.9151 


-30.9388 


18.01 


22.8 


611 


ESQ 434- G 040*^ 


0.6 


Sy2 


R, I, X 


09 47 39.8, -30 56 55 


J0947. 7-1-0725 


146.9447 


7.4191 


5.58 


4.7 


581 


3C 227 


0.4 


Syl 


E.R'^ 




J0959. 5-2251 


149.8805 


-22.8561 


7.80 


10.8 


566 


NGC 3081'= 


1.8 


Sy2 


R 





''Catalogs and Instruments used for identification: R=ROSAT catalog, 1= INTEGRAL catalog, H= HEAO-1 catalog, E= Einstein IPC catalog, X= XRT 
observations. 

''Source identified thanks to its proximity to a ROSAT Bright Source with the method described in the text. 
■^SC 227 was detected in a long (11 ks) ROSAT-PSPC observations jCrawford fc Fabian|[T995h . 

'^Source J081 1.54-0937 has a ROSAT counterpart, IWGA J0811. 5-1-0933, reported in the WGA catalog jWhite et al.lll994l) . This source coincides with RX 
J081132. 4-1-093403. 

'^Object has a ROSAT source in the error box, however the association has to be proven. 
References. — (1) Tueller et al. (2005); (2) Kennea et al. (2005); (3) Ajello et al. (2006); 
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Fig. 1.— Image of LMC X-1 (right) and PSR B0540-69.3 (left) clearly separated. The pixel 
size is 5' and the map is about 1x1 deg^. 



'DEC 
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Fig. 2. — Exposure map, corrected for telescope vignetting, of the survey field presented in 
text after data screening. 
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arcmin 



Fig. 3. — The solid line is the distribution of the offsets of sources detected in the individual 
observations (see Tab. [2]) from their catalog position. The inner histogram (dotted line) 
shows the detections of the Crab while the dashed line is the cumulative distribution of all 
detections. 90% and 95% confidence limits are at radii of 1'.8 and 12 respectively. 
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Fig. 4. — Offset from catalog position of sources detected in individual observations as a 
function of S/N. The solid line is the function described in Eq. [11] 
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Fig. 5. — Significance map of tlie 90x90 deg^ surveyed area. Blue circles mark all sources 
above 4.5 a presented in Tab. HI while red squares show the position of bright sources which 
were included in the background model. 
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Fig. 6. — S/N distribution. The dashed hne is an overlaid Gaussian with a — 1. 




Fig. 7. — Left Panel: Sky coverage as a function of minimum detectable flux for S/N=4.5 in 
the 14-170 keV band. As it can be seen 75% of the surveyed area is covered to 1 mCrab. The 
dashed line is the sky coverage for \b\ > 15°. Right Panel: Pixel 5 a sensitivity threshold 
as a function of pixel exposure time. The dashed line is the fit to data points corresponding 
to the function 3.0mCrab(T/100ks)-° ^ 
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10 100 1000 
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Fig. 8. — Probability of finding at least one ROSAT source within a given radius for our 
source sample of table HI The thick solid line is the probability for the real sample while the 
thin line is the expected contribution of chance coincidences. The inset shows a closeup view 
of the chance coincidence curve for the region of interest (200" < R <300"); the dashed line 
is the chance coincidence distribution for latitudes —40° <b< —20°, while the dotted line is 
for 0° <b< 20°. 




Fig. 9. — Gaussian fit to the curve of Fig. [81 The PSLA at 95% and 99% confidence is 
respectively 3'.3 and 5'. 
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S/N 

Fig. 10. — Offset of the sources detected in the final survey image as a function of their 
significance. The solid hne is the fit to the data with the function OFFSET = 10'.7xS/N~°'^^. 
The long dashed line is the fit to > 3 a deviations from the previous fit and gives the 
maximum expected offset for a given significance. For comparison, the short dashed hne 
shows the best fit (Eq. [11]) to the offset-significance relation for sources detected in individual 
pointings. 



